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ABSTRACT 

Higher-order statistics are a useful and complementary tool for measuring the cluster- 
I ing of galaxies, containing information on the non-gaussian evolution and morphology 

I of large-scale structure in the Universe. In this work we present measurements of the 

three-point correlation function (3PCF) for 187,000 galaxies in the WigglcZ spectro- 
scopic galaxy survey. We explore the WiggleZ 3PCF scale and shape dependence at 
three different epochs z = 0.35, 0.55 and 0.68, the highest redshifts where these mea- 
surements have been made to date. Using N-body simulations to predict the clustering 
of dark matter, we constrain the linear and non-linear bias parameters of WiggleZ 
galaxies with respect to dark matter, and marginalise over them to obtain constraints 
on crg{z), the variance of perturbations on a scale of 8 /i~^Mpc and its evolution with 
redshift. These measurements of asiz), which have 10-20% accuracies, are consistent 
with the predictions of the ACDM concordance cosmology and test this model in a 
new way. 

Key words: cosmology - large scale structure statistics - higher order correlations 



1 INTRODUCTION 



* E-mail: fmarintSastro. swin.edu.au (FAM) 



In t he current str u cture formation p ar adiRm 
(e.g. iPress fc Schechted 1 1974 I White fc Re3 Il978l : 
I White fc FrenkI Il99ll : iBerlind fc Weinberg! |2002| ). galaxies 
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form inside dark matter halos, which evolved from small 
perturbations in the early universe. This allows us to con- 
nect the galaxy field to the overall matter distribution, and 
therefore to use large-scale galaxy clustering to constrain 
cosmological models and their parameters ( s ee for instance 



Peacock etliD I2OOII: lEisenstein et all 120051; IPercival etliD 
20071 : iKazin et al.ll2O10n Sanchez et al.ll2012l and references 



therein). This connection, however, is not a perfect one, 
since galaxy observables such as luminosity, colour, etc. 
are also shaped by baryonic physics and environmental 
effects, with the consequence that diffe rent types of galaxies 
have different clusteriii g properties (jNorberg et al.l I2OOII : 
IZehavi et al.]|2005l . I2OIII ) , described as 'galaxy bias'. 

The galaxy 2-point correlation functions have been the 
main tool to constrain cosmology using large-scale struc- 
ture, because the shape of the 2-point clustering of matter 
depends on cosmological parameters such as the matter den- 
sity, baryon fraction and neutrino mass. In some cases it is 
possible to obtain these constraints marginalizing over the 
bias of the galaxy populations we use; an example is the 
cosmological constraints from Baryon A c oustic Oscillation 
(BAQ) measurements (lEisenstein et al.| 20051: Cole et al 



20051 : iBlake et al.1l2011al : iBeutler et al.ll201ll : ISanchez et al 



2OI2D where parameters such as the cosmic distance scale 



and the Hubble expansion rate H{z) can be measured using 
only the position of the BAO peak, which does not depend 
on the details of the galaxy population used, at first or- 
der. However, there are other cosmological parameters which 
cannot be constrained using 2-point galaxy clustering statis- 
tics only. 

In particular, the amplitude of primordial perturba- 
tions, parameterized in the low-redshift universe by erg (2:), 
the r.m.s. of the matter density field in 8 h^^Mpc spheres ex- 
trapolated to redshift z by linear theory, is degenerate with 
the details of the galaxy population, encoded in the large- 
scale linear galaxy bias parameter bi. On large scales these 
two parameters have the same effect on the overall ampli- 
tude of the galaxy 2-point correlation functions; therefore, 
one can only constrain the product bias. 

Resolving this degeneracy requires use of another ob- 
servable, or adoption of a particular galaxy evolution model. 
For the first approach, other obser vables such as lensing 
(e.g. IFu et al.|[200^ : iLin et aLlbOllI). or the mass fu i iction 
of galaxy clusters (e .g. Eke et alT 19961 : IRozo et al.l I2OI0I : 
iKillpinger eral]|2012l ) have been used to constrain erg; but 
current constraints are degenerate with other parameters 
such as the matter density Qui- The second approach, select- 
ing a particular galaxy population with a known evolution 
of its clustering, allows disentangling of the linear galaxy 
bias and ag, and has been studied in iToieiro et al.l (|2012l ) 
using a passive-evolving luminous subsample of SDSS-I/II 
and SDSS-III survey galaxies. This method gives good con- 
straints on as; but whereas it works well for their galaxy 
sample (of luminous red galaxies), for other galaxy surveys 
it might be difficult to find a suitable passive galaxy popu- 
lation. 

Another, complementary way to break these degenera- 
cies is to measure three-point correlation functions (3PCF) 
of the same galaxy dataset. Two-point statistics are only a 
complete description for Gaussian fields; but the late time 
large-scale structure, driven by no n-linear gravitational clus - 
tering, is strongly non-Gaussian (|Bernardeau et al]|2002h . 



and higher-order correlation functions thus encode addi- 
tional information that can be used to constrain galaxy 
population and cosmological models. The first measure- 
ments of the 3PCF were carried out in angular cat- 
alogues as a way to verify the hierarc hical model of 
structure formation (|Peebles fc Grothlll975l ). and more re- 
cently, in large-scal e spectroscopic survey s such as 2dFGRS 
( Verde et all I2OO2I : Ijing fc Bo rW '2 004 iGaztafiaga et all 
2005h. and SPSS iKavo et al.. .2004|riNichol et al.1 l2006l: 
Kulkarni eral]|2qq7l:TGaztanaga et al.ll2008l : iMcBride et all 



2OIII : lMarfnll201ll : IGuo et aLlbOlSh Trhe main goal of these 



measurements is to test theories of growth of structure and 
the predictions of cosmological simulations, and to measure 
the biasing of the galaxies with respect to the dark matter 
distribution. 

In this work we present results of the measurement 
of the three-point correlation function for a sample of 
187,000 galaxies from the WiggleZ Dark Energy Survey 
(|Drinkwater et a"l]|2010l ) , which probes galaxies in the range 
0.1 < z < 1.0 with a median redshift z --^ 0.6. Using N-body 
simulations to study dark matter statistics, we estimate the 
WiggleZ galaxy bias and thereby measure as- These esti- 
mations have been done in the past jGaztanaga et al.l l2005l : 
IRoss et al.ll2006l : lMan'nll201ll : iMcBride et aLlbOllh . but the 
fact that the WiggleZ survey spans such a large range of 
redshifts with a large overall volume allows us to split our 
galaxy sample in redshift slices and measure as as a func- 
tion of redshift, hence constraining the large-scale structure 
growth history. 

This paper is organized as follows: In §2 we describe 
our survey and the simulations we use; in §3 we introduce 
the galaxy 3PCF and how it is measured, along with the 
model connecting galaxy and dark matter clustering. In §4 
we present the measurements of the WiggleZ 3PCF as a 
function of scale and shapes. In §5 we discuss constraints 
on the galaxy bias and erg as a function of redshift. In §6 
we summarize our findings. We note that a fiducial flat 
ACDM cosmological model with matter density ilm ~ 0.27 
and Hubble parameter Ho = 100 h km s~^ Mpc~^ with 
h = 0.7 is used throughout this paper to convert redshifts 
to distances, which are measured in /i~^Mpc. 



2 DATA AND SIMULATIONS 
2.1 The WiggleZ Galaxy Survey 

The WiggleZ Dark Energy Survey (jDrinkwater et al.|[201ol ) 
is a large-scale galaxy redshift survey p erformed over 276 
nights with the AAOmega spectrograph (|Sharp et al.ll2006h 
on the 3.9m Anglo-Australian Telescope. With a area cov- 
erage of 816 deg^, this survey has mapped 240, 000 bright 
emission-line galaxies over a cosmic volume of approximately 
1 Gpc^ 

Target galaxies in seven different regions were cho- 
sen using UV photomet ric data from the GALEX sur- 
vey (|Martin et al.l l2005l l matched with optical photom- 
etry from the Sloan Digital Sky Survey (SDSS DR4, 
lAdelmam McCarthy etla]|2006l) for regions in the Northern 
Galactic Cap (9-hr, 11-hr and 15- hr), and from the R ed- 
Sequence Cluster Survey 2 (RCS2. ICilbank et al.|[201ll ) for 
those regions in the Southern Galactic Cap (0-hr, 1-hr, 3- 
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Figure 1. Angular distribution of WiggleZ galaxies. The top four regions correspond to those WiggloZ galaxies in the RCS2 footprint; 
the bottom three regions to the ones obtained using SDSS. Colours correspond to sub-regions containing the same effective area. 



Table 1. Number of galaxies in WiggleZ regions used in this 



paper. 


redshift 


[0.1,0.5] 


[0.4,0.8] 


[0.6,1.0] 


JK subregions 


Zeff 


0.35 


0.55 


0.68 




00- hr 


6601 


10698 


8774 


9 


01-hr 


6038 


9437 


7880 


8 


03- hr 


6492 


10241 


8756 


8 


22- hr 


13508 


16146 


11024 


15 


09- hr 


10106 


18978 


11424 


12 


11- hr 


13603 


23940 


13919 


16 


15- hr 


14517 


30015 


19471 


20 


all regions 


74440 


119455 


81248 


88 



hr, 22- hr). The selectio n criteria consisted of ap plying mag- 
nitude and colour cuts (|Drinkwater et a"l]|201(]| ) in order to 
select star-forming galaxies with bright emission lines with 
a redshift distribution centered around z ~ 0.6. The se- 
lected galaxies were observed in 1-hour exposures using the 
AAOmega spectrograph, and their redshifts were estimated 
from strong emission lines. 

To study the evolution of the bias and as with cosmolog- 
ical time, we use three overlapping redshift slices [0.1,0.5], 
[0.4, 0.8] and [0.6, 1.0]. We estimate the effective redshift of 
each sample by averaging the redshifts of galaxy pairs at the 
distances covered by our study, i.e. from 10 to 100 h~^Mpc; 
we find that the effective redshifts for the closest, middle and 
farthest slices are 2eff=[0.35, 0.55, 0.68] respectively. Table 
[T] shows details of the samples used. 

Figure [1] shows the angular distribution of galaxies in 
the regions considered. We show the targets in r.a.-dec. 
coordinates, where it can be seen that the angular com- 
pleteness varies considerably between regions, due to mask- 
ing of bright stars, the availability of input GALEX imag- 
ing, and differences in the accumulated observation time 
within each region and between regions. If not taken into ac- 



count properly by modelling the angular selection function, 
these non- uniformities may lead to artificial structures, dif- 
ferent from what we can expect from cosmic v ariance. Sev- 
eral s t udies, such as [ Ga ztanaga et al.l (r2005f) : Nichol et al.l 
(|2006h : iMcBride et all (|2011); Norberg et all l(201lh . agree 
that higher-order correlation functions are more sensitive to 
these effects than the 2-point function. However, as we de- 
scribe below in §3.2, we conclude that our modelling of the 
angular completeness is adequate to carry out our analyses. 

Figure [5] shows the redshift distribution of the differ- 
ent regions, peaking at z ~ 0.6, but extending to redshift 
z ~ 1.0. The variable number density with redshift deter- 
mines the effective redshifts of our samples as measured 
above. It can also be seen that the average distribution of 
galaxies varies between regions. This is partly explained by 
cosmic variance, but also the selection functions of SDSS and 
RCS2 galaxies differ considerably at low redshifts, owing to 
the avail able colours for galaxy s election from the input cat- 
alogues (jPrinkwater et al.l [2OI0I ') . To deal with these issues 
we model the angular coverage and redshift distrib utions in 
each survey region individually (|Blake et al.ll2010l) . 



2.2 The GiggleZ simulations 

In order to measure galaxy bias we need to model the 
underlying dark matter correlations. For the 2-poi nt func- 
tions there exists a large literature of models (e.g. Peebles! 



198C 



2002 



Kaiseij|l987l. ICoorav fc Shethll200i iBernardeau et al.l 
Smfth et alll2003l , and references therein), but in the 



case of the higher-order correlations, modelling has focused 
mostly on the large-scal e behaviour fe.g ljing fc Borne^ll997l : 



iBernardeau et al.|[20o3 ). although there have been efforts 
to model the non-l inear, small scales (e.g lYang et al.|[2002l : 
iFosalba et"al]|2005^ ■ Most importantly, for the 3PCF there 
is no satisfactory treatment of redshift-space distortions, al- 
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Figure 2. Rcdshift distribution of WiggleZ galaxies. Thin linos 
correspond to the radial selection function for each individual 
angular region. Thick lines represent the redshift distribution of 
all WiggleZ galaxies (black), WiggleZ galaxies in the North Cap 
(blue) and WiggleZ galaxies in the South Cap (red). 



though some attempts have been made f or the bispectrum, 
the F ou rier transform of th e 3PCF, by IScoccimarro et ahl 
l| 19991 ): ISmith et all (|2008l '). which are vaUd for limited 
range of scales. There f ore, a s has been d o ne in previous 
workslGaztanaga et al.l (|2005l ): [Marfnl (|201ll ): lMc"Bride et all 
l|201ll ). we will obtain constraints on the galaxy bias and 
cosmological parameters by comparing the WiggleZ galaxy 
3PCF with the dark matter correlations measured in N-body 
simulations, which include the full set of non-linear effects. 

We measured the dark matter correlation functions us- 
ing the Gigapars ec WiggleZ Survey simulations (GiggleZ, 
IPoole et al.ll2012l ). which have been generated in support of 
WiggleZ science. In a 1 /i~^Gpc^ periodic cube, 2160^ dark 
matter particles with individual masses of nip = 7.5 x lO'' 
were evolved using a fiat ACDM r nodel, with cos- 
molo gical parameters from WMAP5 results (jKomatsu et al.l 
I2OO9I ). In order to compare the correlation functions of Wig- 
gleZ galaxies and dark matter, we measured them in snap- 
shots of the simulation at the same effective redshifts as the 
WiggleZ subsamples. 

Dark matter halos are identified in two steps 
(jSpringel et al.|[200ll ): Firstly, using a friends-of- friends algo- 
rithm with a linking length of I — 0.2 times the mean particle 
separation, bound structures are found (parent halos). Sec- 
ondly, given the high resolution of our simulation, we were 
able to find gravitationally bound substructures inside these 
parent halos. From the main halo catalog, we create subsets 
ordered by maximum circular velocity (which we use as a 
proxy for halo mass) with the same mean number density 
as the WiggleZ galaxies n ~ 2.5 x 10~* (/i^^Mpc)"^. These 
halo catalogues are used to carry out consistency checks in 
our phenomenological model to estimate the bias and cos- 
mological parameters. 



3 THE GALAXY THREE-POINT 
CORRELATION FUNCTION 

3.1 Definitions and Methods 

The galaxy n-point correlation functions are the average of 
correlat ed galaxy ove rdensity 5gai measured at n different 
points (|Peebleslfl980l ). Whereas the two-point correlation 
function ^(r) (2PCF) allows us to estimate the probability 
of finding pairs with a separation r, the three-point correla- 
tion function (^(ri, r2, rs) (3PCF) describes the probabilities 
of finding triplets with galaxies as vertices. The joint proba- 
bility of finding three objects in three infinitesimal volumes 
dVi, dV2, and dV-j is given by the 'fuU' three-point correla- 
tion function: (jPeebles fc Grothll 19751 : |Peebleslll980l ) 



[1 + ^{ri) + C(r-2) + ^{r:i) + C(ri, ra, rg)] 



xn^dVidV2dV3 



(1) 



where n is the mean density of objects, ^ is the 2PCF, and 
is the reduced or 'connected' 3PCF . In other words, this 
means that the probability to find galaxies in a particu- 
lar triangular configuration has contributions from triplets 
found by random chance, plus contributions from correlated 
pairs plus the third point found at random (the ^ terms), 
and lastly by intrinsically correlated triplets (the ^ term). 
In the dark matter or galaxy field the 2PCF and 3PCF are 
given by: 



C{ri,r2,ri) 



= (5(x05(x + r)) 
= (5(xi)5(x2)5(x3)), 



(2) 
(3) 



where S is the fractional overdensity of objects (galax- 
ies, halos or dark matter particles) or the continuous 
field studied, and xi, X2 and xa form a closed trian- 
gle (see Figure |3]). The triangle sides are the dis- 
tances between objects in the triplet; thus the 3PCF de- 
pends upon the scales and shapes of spatial structures 
dBarriga fc Gaztafiaga 20021 : Sefusatti fc Scoccimarrall2005l : 



iGaztafiaga fc Scoccimarrd l2005l : iMan'n et al.1 l2008h Since 
the ratio C/?'^ is both predicted on larg e scales from per- 
turbation theory (jBernardeau et aLll2002l ) and measured to 
be close to unity over a large range of length scal es, even 
thoug h ^ and each vary by orders of magnitude (jPeeblesI 
ll98Ci ). we will often present results using the 'reduced' (or 
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catalog o f the region, respe c tively . The 3PCF is calculated 
using the lSzapudi fc Szalavl (|l998l ) estimator: 




SyS^ = 20:40 h ^Mpc 




Figure 4. Correlation function measurements for a GiggleZ halo 
catalog. Selecting a set of triangular configurations, with s = 15 
/i~^Mpc (top) and s = 20 h^^Mpc (bottom) and u = 2, we 
plot the redshift-space reduced 2PCF of the third side i{s3{6)) 
(left panels) and connected 3PCF, (right panels) of a selected 
DM halo catalogues from the GiggleZ simulation with similar 
clustering as WiggleZ galaxies, showing the effects of radial and 
angular selection functions. 



normalized) 3PCF, 
Q{s,u,0) = 



C(s,u, e) 



Here, s = ri sets the scale size of the triangle, and the 
shape parameters are given by the ratio of two sides of the 
triangle, u = r2/ri, and the angle between those two sides, 
9 = cos~^(fi • f2), where fi, r2 are unit vectors in the direc- 
tions of those sides. The reduced 3PCF is also better suited 
for visualizing the growth of non-gaussian structure and the 
shape dependence of clustering than On the other hand on 
large scales, since ^ — >■ 0, the ratio C/C^ becomes very un- 
stable and its errors non-Gaussian, with the consequences 
of overestimating the covariances, diminishing the overall 
signal-to-noise ratio and introducing a systematic deviation 
in the confidence intervals of the fitted parameters, therefore 
we use ( and not Q for model fits. 

In this work we measure correlation functions for trian- 
gles with base side s = 10, 15, 20 and 30 h^^Mpc, with the 
shape parameters u = 1.0, 2.0, 3.0 and 10 equally spaced 
bins in 9. In total we measure correlations for 120 triangular 
configurations. 

3.2 Measuring Correlation Functions 

We measure first the 2PCF and 3PCF in each WiggleZ re- 
gion (i.e in angular and redshift cuts). For a particular Wig- 
gleZ region, we c alcula te the 2PCF using the estimator of 
iLandv fc Szalavl (|l993h . 



DP - 2DR + RR 



(4) 



Here, DD is the number of data pairs normalized by Nd x 
No/2, DR is the number of pairs using data and random 
catalogues normalized by NdNh, and RR is the number of 
random pairs normalized by Nr x Nit/2, where No and Nr 
are the number of points in the data and in the random 



ODD - 3DDR + 3DRR - RRR 

Rim 



(5) 



where DDD, the number of data triplets, is normalized 
by No/6, and RRR, the random data triplets, is nor- 
malized by iV^/6. DDR is normahzed by NoNr/2, and 
DRR by NdNr/2. Due to the low density of our galax- 
ies n ~ 2.5 X 10"" h/^ Mpc"^, we are limited by shot noise 
on these large scales and consequen tly the application of 
FKP weighting (|Feldman et al.lll994 ) to the pair and triplet 
counts does not affect the results. To estimate the number of 
pairs and triplets, we use the ntropy-npoint software, an ex- 
act n-point calculator which uses a kd-tree framework with 
true parallel capabil i ty and enhanced rou tine performance 
^Gardner et al.ll2007l : [McBride et al.ll201ll ). 

The ran dom catalog u es we re built using the methods 
described by iBlake etlo] (|2010l ), which estimate the angu- 
lar and radial selection function of each survey region due 
to survey geometry and incompleteness in the parent pho- 
tometry and spectroscopic follow-up. This modelling process 
produces a series of Monte Carlo random realisations of the 
angular and redshift catalogue in each region, which is used 
in our correlation function estimations. In this paper, we 
measure the 3PCF using 10 random catalogues for each re- 
gion, with Nh — 4Nd for each of the random catalogues for 
the intermediate scales (s = 10, 15 /i~^Mpc), and Nr = No 
for the largest scales (s = 20, 30 h~^Mpc). For our choice of 
bi nning (resolut i on) of the triang l es, we use the same scheme 
as iMarinl (|2011^ : [McBride et all (120111 ): first, we select the 
central s, u and 9, and their corresponding side lengths in 
redshift space Si, with i — 1,2,3. Then, to calculate the 
2PCF and 3PCF, we accept triangles with sides between 
(1 - 0.075)s, < n < {1 + 0.075)si, implying a 15% bin- 
ning re solution. Thi s is a higher resolution than used for the 
LRGs (|Man'nll201]] ) but is justified by the higher number 
density of the WiggleZ galaxies. 

In FigureUwe explore the effects of the radial and angu- 
lar selection functions on the 2PCF and connected 3PCF ^ 
for a selection of configurations (s = 15, 20 /i~^Mpc, u = 2). 
Using a dark matter halo catalog from the GiggleZ simu- 
lation that has a similar 2-point clustering as our WiggleZ 
sample ai z = 0.55, we create three different kinds of mock 
catalogues with the same geometry as the survey: the first 
group (in red, dashed line) does not include any radial or 
angular selection effects apart from the large-scale bound- 
aries of the WiggleZ regions. Green (dotted-dashed) lines 
are measurements from mocks with the same radial selec- 
tion function of WiggleZ galaxies at ZeB = 0.55. Blue (solid) 
lines are mocks with the same angular and radial selection 
function of WiggleZ galaxies. In general, the measurements 
of the correlation functions using the different mocks are 
very similar, signaling that our selection functions and ran- 
dom catalogues allow us to recover the intrinsic correlations 
of the galaxy field. 



3.3 Galaxy Bias Model 

Since different types of galaxies form inside different dark 
matter halos, they are an imperfect tracer o f the overall dark 
matter distribution (jBardeen et al.l Il986l : ICoorav fc ShethI 
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\2OOi : iBerlind fc Weinberd [200^ '). and their n-point cor- 
relations will differ as well. Many models of this galaxy 
bias have been proposed, and an accepted working model 
on large scales is the deterministic and local bias formal- 
ism, where we relate real space ga laxy overdensity ^gai to 
the underlying matter dens ity Sm ()Frv fc Gaztanagalll993l : 
iFrieman fc Gaztanagalll994 ): 



real space 



(6) 



up to second order, where Sg^i and 5m are the local galaxy 
and matter overdensities smoothed over some scale R. To 
leading order, this bias prescription leads to a relation be- 
tween th e galaxy and matter 2 PCF and connected 3PCF. 
Following TajT&Szapud? 1*2005') , leading order Perturbation 
Theory (|Bcrnardcau ct al...200Z ) shows that if we fix all cos- 
mological parameters except the overall amplitude of the 
initial spectrum of perturbations characterized by as, then 
there is a degeneracy between the effect of this parameter 
and the bias on the 2-point and 3-point functions. The rela- 
tions between matter and galaxy correlations in this model 
axe: 

Cgal(7') = (o-8/o-8,fid)^fo?Cm(r) (7) 

62fe?(Cm(?'-i2)5m(r23) -|-perm.)]. (8) 

In observations we measure the correlation functions in red- 
shift space, where the 'real' correlations are distorted by 
peculiar velocities (with respect to the Hubble expansion), 
which on large scales depend on the growth rate of pertur- 
bations / « Qm{z)^'^^ and on galaxy bias. For the 2-point 
funct ion we have that ^z-spaco ^ /2^r-spaco, with (jKaised 
1 19871 ): 



(9) 



In the case of the 3PCF, there is also an effect from red- 
shift space distortions of similar order jScoccimarro et al] 
1 19991 ). However, it depends not only on the linear bias and 
/, but also on the non-linear bias and the shape and scale 
parameters of the triangle observed. As mentioned before, 
analytical models of these distorti ons have been proposed 
for the bispectrum on large scales (|Scoccimarro et al.lll999l ) 
where its validity is limited. On small, non-linear scales 
(ISmith et al ]|2008l ) the transformation back to configuration 
space is challenging (six-dimensional integrals) and numeri- 
cally intractable. 

For these reas ons we opt to us e an empirical model that 
has b een used by |Pan fc Szapudil (|2005h ; iGaztafiaga et all 
(|2005l ) in analyses of the 2dFGRS galaxies. We use the N- 
body measurements of the correlation functions in redshift 
space in equations (7) and (8), replacing ^r-spacc —>■ Cz-spacc 
and Cr-spacc — ^ Cz-space at the different effective redshifts 
of our WiggleZ slices. Given the low signal-to- noise ratio of 
the 3PCF measurements, this is justified by the fact that for 
low bias tracers (such as the 2dFGRS galaxies studied in the 
papers mentioned above) with 6 ~ 1, f2^bi ~ /2,6i=i, with 
differences of the order of 10% when (fei — 1) = 0.3, smaller 
than the statistical error in our measurements (of the order 
20%-30%); the impact on the constraints on as is slightly 
lower. Since we are using N-body simulation measurements 
to compare to our galaxy correlations, the only cosmological 
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Figure 5. The 2PCF ?(s3(6»)) (loft) and connected 3PCF ( (right) 
of two samples (groups) of halos from the GiggleZ simulations 
in real space (top) and redshift space (bottom) at a snapshot 
z = 0.6. Black lines are results from the dark matter particles 
from the simulation, blue triangles display measurements for a 
halo sample ('group 05') with average mass 4.8 X 10^^ h~^MQ; red 
squares show results for a halo sample ('group 44') with average 
mass 6.9x 10^^ h~^MQ. Blue (short dashed) and red (dot dashed) 
arc biased DM for low mass and high mass halos respectively (bias 
parameters are listed in the figures). 



parameter we modify is as through the ratio as/agfid in eqs. 
(7) and (8). 

The justification of this model is illustrated by Figure 
[S] where we take two halo samples (groups) from the Gig- 
gleZ simulation at redshift z — 0.6, with a similar number 
density as WiggleZ galaxies Uh = 2.5 x 10^'' (ft^^Mpc)"^. 
The first sample, 'group 05', is composed of low mass ha- 
los, with clustering similar to that of WiggleZ galaxies, and 
the second sample, 'group 44', consists of very massive ha- 
los. We measured their correlation functions in real space 
and then estimated their linear and non-linear bias param- 
eters by comparison with the matter correlation functions 
measured in the GiggleZ simulation, using eqs. (7) and (8) 
for a fixed as{z = 0) = 0.812 (see §5.1). In the case of the 
low bias sample, Xr.g05/dof= 0.87, and for the high mass 
sample Xr,g44/dof= 0.96 We then performed measurements 
in redshift space; we observe that in the case of the low- 
bias tracers, using the same bias in redshift space fits well 
their redshift space 2PCF and 3PCF, with xl,g05/'ioi= 1.07, 
but the same can not be said for the high-mass halos where 
Xz,g44/dof= 22.44. Since WiggleZ galaxies have a low l inear 
bias 6i ~ 1 (|Blake et al ] |2011bl : [Contreras et al]|2012l ). the 
approach we take is adequate for obtaining measurements of 
the linear and nonlinear bias parameters and erg as a func- 
tion of redshift. 



4 RESULTS 

In this section we present the measurements of the 2PCF ^, 
connected 3PCF ({s,u,9) and reduced 3PCF, Q{s,u,6), of 
WiggleZ galaxies for a range of scales and shapes at different 
redshifts. We explore differences between regions, evolution 
with redshift, and our estimation of statistical errors and 
covariance. 
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Figure 6. The 2PCF ^(sa) (top left), the connected 3PCF 
C (middle left) and the reduced 3PCF Q (bottom left) of 
each WiggleZ region in the Zeg = 0.55 slice, for triangles 
with SI = 15 h~^Mpc and ti = 2. In the right-hand plots 
we show the corresponding diagonal errors. 

4.1 Building the WiggleZ Survey correlations 

In Figure [6] we show measurements of the correlation func- 
tions in the ZeS ~ 0.55 sample for each WiggleZ region. Mea- 
surements in different regions are consistent within the sta- 
tistical errors from cosmic variance and shot noise. On these 
scales we notice how small differences in ^ and ( translate to 
larger discrepancies in Q. The noisiest 3-point functions are 
obtained in the smallest regions (in terms of volume), in this 
case the regions overlapping with the RCS2 survey. We will 
build a 'combined' set of correlation functions calculated by 
optimally weighting individual contributions of the regions. 

To measure the diagonal errors and covariance matri- 
ces, ideally we require a large number of mock galaxy cat- 
alogues whose correlation functions have similar amplitude 
and shape dependence to the one observed in our data. How- 
ever, this is not currently tractable for WiggleZ galaxies 
due to their low bias, which would necessitate many high- 
resolution simulations in a cosmological volu me. In past 
studies we have used lognormal realisations (|Blake et al.l 
to generate a large number of mocks suited to match 
the WiggleZ 2-point amplitudes. Unfortunately, by their 
construction (from generating points with a particular 2- 
point distribution), these are not capable of reproducing the 
higher-order clustering of galaxies. 

In this work, error measurements in each region are 
calculated from jack-knife resampling (|Zehavi et al.l l2005l : 
iMcBride et al.ll201ll : lNorberg et al ]|201ll ). In this method we 
divide the whole volume of the sample in identical subsam- 
ples i — 1... A'^, and we then measure correlation functions 
for the whole volume minus the i-th subsample times to 
get a set of A'' correlated measurements. In our case, in all 
WiggleZ regions we take subregions of equal area (weighted 
by sky completeness), with an equivalent physical size of 
approx. 120 x 150 x 900 {h~^Mpcf a.t z = 0.55. The an- 



gular cut has the same area independent of region, as can 
be seen in Figure [1] thus some WiggleZ regions have more 
JK subregions than others, depending on their total area 
coverage. 

In order to obtain the JK vaxiance in each region, we 
measure each Xi statistic (where Xi can be the 2PCF or the 
3PCF), subtracting one of the JK subregions in turn. Then 
we calculate the variance axi of the individual WiggleZ re- 
gion as 



2 



TV/ 



N: 



5^ ^(Xf,,-(X,>^); 



(10) 



where N^'^ is the number of jack-knife subsamples, in re- 
gion i (see Table 1). Then, we calculate the correlations in 
the overall survey using inverse- variance weighting. For the 
statistic Xcomb, this is calculated as 
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E 



Xi 



(11) 



where ax is the variance of the statistic in the WiggleZ sub- 
region (calculated in equation 10), taken from the jack-knife 
resampling method. Nreg = 7 is the number of WiggleZ re- 
gions we use for the calculation of the combined correlations. 
We do this for ^, ( and Q. In principle there should be no dif- 
ference between calculating correlations using this method 
and measuring triplet counts across the whole survey; in 
practice our method is more computationally efficient and 
gives us extra systematics tests by allowing us to compare 
results region by region. 

Overall diagonal errors and covariance matrices are cal- 
culated by jack-knife resampling the whole set of survey re- 
gions (see Figure [1]). This way, we produce a catalogue of 
NtJt = 88 measurements. The variance of the correlations 
is calculated as 
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N. 



ATJh 



-E(^: 



{XY 



(12) 



3=1 



We also use this method to calculate the covariance matrix, 
which is used in the Maximum Likelihood approach to mea- 
sure the galaxy bias and cosmological parameters. 

4.2 The Combined WiggleZ 3PCF 

Figure [7] shows the measurements of the redshift-space 
2PCF C(s3(s,«,e)), connected 3PCF C(si,S2,6') and re- 
duced 3PCF, Q{si, S2,d), of WiggleZ galaxies (from opti- 
mally combining the seven independent regions) in all red- 
shift slices for a range of scales si = 10, 15, 20 and 30 
/i~^Mpc and shape u = 2 as a function of 9. We have addi- 
tionally measured the correlations using u = 1 and = 3, for 
a total of 120 configurations, which are shown in Appendix 
A. 

We note qualitatively that we recover the same shape 
dependence of the galaxy 3PCF (mostly noticeable when 
looking at Q) which has been observed for galaxies at low 
redshift (|McBride et al.ll201ll : iMan'nl l201ll '). namely a big- 
ger 3PCF amplitude at small and large 9, i.e the collapsed 
and elongated configurations. This 'V-shape' is more promi- 
nent for large scales and elongated shapes; it is a conse- 
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Figure 7. The combined redsliift-space 2PCF C(s3(6»)) (left), connected 3PCF C,{s,u = 2, 6) (middle), and reduced 3PCF Q{s,u = 2, 0) 
(right) of WiggleZ galaxies in the Zgfj = 0.55 slice (blue triangles and error bars), in the slice at z^g = 0.35 (green dashed line) and in 
the Zgff = 0.68 slice (red dashed-dotted line). Different rows cover the range of scales of the triangles. Errors have been determined by 
jack-knife resampling. 



quence of the morphology of galaxy structures varying from 
spherically-shaped clusters and groups on small ~ 1 h^^Mpc 
scales to filaments on the largest scales. This shape depen- 
dence of the 3PCF depends on the galaxy type under inves- 
tigation. It has been observed in S PSS (jKavo et aD l200 
iMcBride et~alll201ll ) and 2dFGRS (jCaztanaga et al-ilm 
that L* blue galaxies tend to have small 3PCF amplitudes on 
small scales and very pronounced V-shapes on large scales, 
compa red to red g alaxies and to L > L* galaxies, such as 
LRGs l|Marfnll201ll ). Although at a different redshift, we find 
that the shape dependence of the WiggleZ 3PCF agrees with 
these lower redshift measurements. 

We also note that on larger scales the behaviour of the 
reduced 3PCF Q{6) for the most elongated shapes is more 
erratic, specially for the z — 0.35 and z = 0.68 redshift 
slices. This is due to the fact that as ^ — > on large scales, 
the measurements of Q are less robust and its errors become 
non-Gaussian. But if we turn to analyze ( instead, we can 
see that it is adequately measured up to the largest scales 
shown here, sg ~ 100 h'^Mpc. 

Comparing the clustering signal in different redshift 
slices, we can see that in general the differences in the 3PCF 
are small and the signal is weaker than in the case of the cen- 
tral Zcft = 0.55 slice. This does not necessarily indicate that 
there is no evolution of the clustering of WiggleZ galaxies 
with redshift; the underlying dark matter clustering changes 



with redshift, and consequently the linear and non-linear 
bias factors evolve. From the 2PCF measurements we can 
estimate the evolution of the linear bias, and using the 3PCF 
we can also test if there is evolution in the non-linear bias 
parameter. 

In Figure [8] we illustrate how the errors in our mea- 
surements vary with redshift by showing the l-cr diagonal 
errors (from jack-knife resampling) of the 2PCF and 3PCF 
measurements for a selection of configurations (s = 15, 20 
/i~^Mpc and = 2) of our (combined) redshift samples. It 
can be clearly seen that measurements are more accurate 
in the central redshift slice (zcH — 0.55) than in the outer 
ones. In these configurations, the relative error in the 2PCF 
is around f^/^ ~ 0.15, for the connected 3PCF the relative 
errors reach o"c/C ~ 0-5 in the central redshift slice. 



4.3 Covariance estimation 

We estimate the correlations between measurements of the 
2PCF 5 and 3PCF ( by empirically calculating the covari- 
ance matrix. Using the jack-knife method, given a number 
of measurements Njk in number of bins Nb, a fractional 
error of a quantity X for the sample k can be written as 



A? 



(13) 
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Figure 8. Diagonal errors of the correlation functions from jack- 
knife resampling. Top: errors for the s = 15 /i^^Mpc, u = 2 
triangles for the 2PCF tr^ (left), connected 3PCF a^^ (centre), and 
reduced 3PCF aq (right); solid line corresponds to measurements 
in the z^ff = 0.55 slice, dashed line for the z = 0.35 slice, and 
dotted-dashed line for the z = 0.68 redshift slice. Bottom: same 
quantities for the s = 20 /i~^Mpc, u = 2 triangles. 
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Figure 9. The correlation matrix of the WiggleZ 2PCF and 
3PCF of triangles with si = 15 h-^Mpc with it = 1, 2 and 3. 
Each element of the matrix is the covariance of each s , u and 9 
triplet. 



where in our case, if i ^ Nt then X = ^(sa) and otherwise 
X = ({si, S2, 83(8)), and axi is the standard error on Xi 
calculated using the jack-knife method. Then we calculate 
the correlation matrix (covariance matrix normalized by di- 
agonal errors) as 



Cij — 



Nj 



k \k 
^3 



(14) 



As an example that shows the observed behaviour of all con- 
figurations at different redshift samples, we show in Figure 
[9] the correlation matrix of both ii^s^iO)) and C,[6) for the 
configurations with s = 15 /i~^Mpc, and it = 1, 2 and 3. In 
this case the number of bins Nb ~ 30, making dj a 60 x 60 
matrix. We divide this matrix into four regions depicting the 
auto- and cross- correlations. In the case of the ^(53) — ^(^3) 
covariance we notice that although it is dominated by di- 
agonal terms, the off-diagonal terms are important too; the 
white stripes (signaling high covariance) in the off-diagonal 
matrices (^(^1,83) — ^(112,53) where ui / U2) are triangles 
that share a similar S3. The ("-correlation matrix also shows 



important non-diagonal elements that are more correlated 
for the elongated shapes (a combination of true covariance 
and binning). The ^(53) — cross-covariance is small but 
needs to be considered in the analysis. 



5 CONSTRAINTS ON GALAXY BIAS AND as 

We compare the dark matter correlations measured in the 
GiggleZ simulations to the WiggleZ 2PCF and 3PCF in or- 
der to constrain the linear and non-linear bias parameters 
using the local bias model described in §3.3. In this anal- 
ysis we assume all cosmological parameters are fixed, and 
fit for as{z) by scaling the amplitude of the dark matter 
correlations in the manner of equations (7) and (8), with all 
quantities measured in redshift space. We want to emphasize 
that since we are using N-body simulation measurements 
to compare to our galaxy correlations, and not an analyt- 
ical model, the only cosmological parameter we can mod- 
ify is as{z) through the ratio crg/o-g^fid, where the fiducial 
value corresponds to the one used in the GiggleZ simulation 
as,fid(2 = 0) = 0.812. 



5.1 Methods 

We carry out a Maximum Likelihood parameter estimation 
test, where we look to minimize the quantity 



X 



i=2JVi, j=2iVb 

E E 

i=i j=i 



(15) 



where A''^ is the number of triangular configurations used; 
we have Nt distances where we measure the 2PGF of S3, and 
Nb triangles where we measure the 3PCF, therefore we have 
(2Nb)^ elements in our covariance matrix. The value of Ai 
is the difference between the correlation measured and the 
biased DM correlation: 



C(s,n,e)r'''=')KW' iori>Nil7) 



A, = (C(s3)f ^ - C(s3)r°*=')M{.3)». for i < Nb (16) 
A, = (C(s,it,e)° 

where C"'*''' and C™"'*''' are given by eqs. (7,8). 

In or der to invert the covariance matr ix we use the ap- 
proach of lGaztafiaga fc Scoccimarrd (|2005|) and repeated in 
several 3PCF work s ( Gaztaiiaga et al. 20051 ; iMcBride et al.l 
I2OIII ; lMari'nll201l] ). of employing only the highest eigen- 
modes of the covariance matrix to minimize effects of nu- 
merical noise. We employ the Singular Value Decomposi- 
tion method, where our normalized covariance matrix can 
be decomposed C = UDV^ (and where V — U for a 
symmetric matrix), where the diagonal matrix Dij = XiSij 
stores the eigenvalues in decreasing order, and the columns 
of the matrix U stores the eigenmodes of C . When invert- 
ing the matrix C"^ = FD'^C/"^, (where D".^ = (1/Af)5„) 
we need to discard some of these eigenmodes, meaning set 
some D~i^ =0. Firstly, we use a finite number of jack-knife 
samples to estimate our covariance. Since using jack-knife 
samples assumes the statistical independence of the sub- 
samples, the jack-knife regions should have a larger spa- 
tial extent than the largest scales studied. In our case we 
use NjK = 88 JK regions because our largest scales are of 
the order 100 /i~^Mpc. However, we are using a large 
number of 2PCF and 3PCF measurements, and generally 
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Figure 10. Dependence of the best-fitting bias parameters and 
trg on tile number of eigenmodes used for tlie z^g = 0.55 2PCF 
and 3PCF analysis. In tlie left panel we show eigenmodes used af- 
ter SVD (from the > y/2/NjK limit), with fixed S3,min = 20 
and ss^max = 120 h~^Mpc. The top left panel shows the total 
S/N according to equation l|18|l as a function of number of eigen- 
modes used. The top right panel shows the x'^/dof, and lower 
panels show the 68% CL intervals for bi, 62 and erg. The vertical 
dashed line represents the eigenmode limit where A? = ^ 2 /Njjf , 
with NjK = 88. 




Figure 11. Dependence of the best-fitting bias parameters and 
(Tg on the fitting range S3 mini with ss^max = 120 fc-lMpcfor the 
Zgff = 0.55 2PCF and 3PCF analysis. The top left panel shows the 
number of eigenmodes used where A? = ^ 2 /Njk as a function 
of S3 mill- The top middle panel shows the total S/N, top right 
panel the x^/dof of the best fit parameters. Lower panels show 
the 68% CL intervals for bi, 62 and erg. 

NjK < Nb, which will make our matrix singular for modes 
i > NjK — 1 (see P ress et al.i il992 for instance). There- 
fore, all those eigenmodes have to be discarded. A second 
cut comes from the fact that even though the covariance 
matrix is not singular when the first cut is applied, there 
are still eigenvalues Af with very low numerical value, which 
will make unstable the inversion of C^J with A? < JK 
iGaztafiaga fc Scoccimarrcl |2005| ). Adding these unstable 
modes, as explained by these authors and later in section 
5.2, is equivalent to introducing artificial 'signal' to our mea- 
surements that will bias our fits. 
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Figure 12. Same as Figure 11, in this case fixing S3,min=20 
/i~^Mpc and having quantities shown as a function of S3, max- 

We also have to set the minimum and maximum scale of 
the model fit. In our analysis this means that we chose con- 
figurations with the third size in a range S3 — [s3,min, -53, max]. 
The minimum scale is given by the range of validity of the lo- 
cal bias model. The maximum scale could be set by system- 
atics in the selection function or when the correlation signal 
is weak. In Figures IIOI 1111 and [12] we show how our choices 
of the number of eigenvalues used, the values of S3, mill and 
S3, max respectively affect our constraints for our z^g = 0.55 
sample. 

We wish to make a sensible default choice for these op- 
tions and establish that our essential conclusions are not 
very sensitive to this choice. Figure [10] investigates the 
dependence of our results on the number of eigenmodes 
included, fixing S3,miii ~ 20 h^^Mpc and S3, max = 120 
/i~^Mpc. We observe that if we use a small number of eigen- 
values we have large variations in the best fit values and poor 
constraints, especially in 62 ; then there is a range when the 
number of eigenvalues used is ~ 50 where our results are in- 
sensitive, and when we include modes with > \pTfWjK ~ 
65 we again have unstable behaviour. In this figure we also 
consider the dependence of the minimum ^ per degree of 
freedom on the number of eigenmodes employed, where the 
degrees of freedom are equal to the eigenmodes used minus 
the number of parameters we seek to constrain. We note 
that good fits are produced for a wide range of choices. We 
can also estima te the total signal-to-noise (S/N) ratio of the 
modes used as (iGaztafiaga fc Scoccimarroll2005h : 




where Xj = (^j when j < Nt, and Xj = (j~N^ when j > Nb- 
Our conclusion from this analysis is that any systematic 
fluctuations in our parameter measurements do not domi- 
nate the statistical errors for a wide range of choices. Our de- 
fault fits are performed for an eigenmode cut A^ — ^/2/Njk 
and a fitting range S3,min = 20 /i~^Mpc and S3, max = 120 
h~^Mpc for the z^g = 0.55 and 0.68 slices, and S3, max ~ 100 
h~^Mpc for the ZcS ~ 0.35 sample. 

Another aspect to consider in our analysis is that the 
measured redshifts in our survey cont ain a small fraction of 
'redshift blunders' (|Blake et al.il2010l ). failures that tend to 
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Table 2. Constraints on bias parameters and erg for WiggleZ 
samples 



Zeff 


6i 


b2 






S/N 


0.35 
0.55 
0.68 


QO+0.10 

1 06+°-^^ 


-0-36Io.o8 
-0.48t°:l* 


0.691-.- 
0.53t°:«? 


1.10 
0.96 
0.82 


3.25 
4.99 
4.62 



wash out the clustering we measure in our galaxy field. This 
redshift blunder fraction is /(, = 0.03 for the Zeff = 0.35 and 
0.55 redshift slices, and = 0.05 for the outer z^s = 0.68 
slice; the correction to the 2PCF (since ^ oc (5gai5gai)) is 
to multiply the data and errors by (1 — /j,)^^, meaning that 
the clustering amplitude increases; and the correction to the 
3-pt function (where C, oc (iJgai'Jgai'Jgai)) is (1 — fb)^^- 




5.2 Constraints at z^g = 0.55 

In Figure [13] we show the measurements of 6i, 62 and ag in 
the Zctf = 0.55 redshift slice. The measured linear bias, 61 ~ 
1, agrees with values obtai ned for WiggleZ ga laxies using 
other methods (|Blakc ct aLll2011bl : IContreras et al. 201^. 
using 2-point statistics, where in our study we additionally 
marginalise over all values of as- We detect a significantly 
non-zero value for the non-linear bias 62 ~ —0.4. 

We measure the amplitude of fiuctuations erg (2) with 
10% accuracy and find that our results agree with indepen- 
dent predictions, based on cosmological parameter measure- 
ments from the CMB in a ACDM model. It is important 
to note that these estimates are independent of any other 
observable than the galaxy clustering itself. Extrapolated 
to z = 0, our measurements of erg from WiggleZ galaxies 
are consistent with conclusions fr om 3PCF m easurements 
of other tracers such as the LRGs (|Man'nll20lil ). which pro- 
vides evidence that this method is robust against the type 
of galaxy used. 

As is shown in Table [2l the empirical bias model is an 
adequate fit to the data, and that can be graphically seen in 
the right-hand plots of Figure [T31 where we show the biased 
dark matter Q(6), which in our model depends on the bias 
parameters but not on ag: 



1 



+ 



(19) 



We can see that our galaxy bias model is adequate on all 
scales, but our fits are mainly driven by the s = 10 and 
15 h~^yipc configurations, which have the highest signal- 
to-noise ratio. 



5.3 Constraints at different redshifts 

We repeated the analysis of the correlation functions for the 
other two redshift slices in order to get constraints on the 
bias parameters and erg as a function of redshift. In Table 
[2] we show our results. In general we can see that there is 
clear trend in all the parameters with redshift, and that 
from the values of x^/dof, our model of the bias provides 
a good fit. We also estimated the covariance between the 
best-fit parameters in different redshift slices, expected due 
to the overlap between redshift samples. In order to make 
this measurement we fitted the bias parameters and ag to 
each delete- 1 jack-knife sample, obtaining Njk = 88 sets of 



Figure 14. Correlation matrix of WiggleZ best fit parameters bi , 
62 and (Tg as function of redshift 



best fit parameters which we used to construct a covariance 
matrix. We show the correlation matrix of these parameters 
in Figure [141 We observe that in any individual redshift slice 
there is a positive correlation between 61 and 62 and a nega- 
tive correlation between these parameters and ag. Although 
there is overlap in the redshifts of the samples studied, the 
correlation coefficients between best-fit parameters in sepa- 
rate redshift slices are small. 

In the following we study the change of the bias param- 
eters and erg (2:) with redshift in more detail. 



5.3.1 Constraints of WiggleZ bias 

Figure [15] displays the change of the galaxy bias with red- 
shift. Due to our selection criteria, we are not necessarily 
selecti ng the same p opulation of galaxies at different red- 
shifts (|Li et al.ll2012l ). specially at redshifts z < 0.5. There- 
fore we only can make general statements about evolution 
of the bias of WiggleZ galaxies. For the linear bias fei it can 
be seen that there is an upward trend of bias with redshift, 
consistent with what is expected for the evolution of halos 
of a given mass: massive objects are less common in older 
times, and therefore more biased. 

In order to compare our measured bias evolution to the 
expectation of simple models, we also plot the evolution of 
the bias of dark matter halos of fixed mass with redshift 
in the halo model. These are given by the a nalytical ex- 
pecta tion using Sheth-Tormen mass functions (jSheth et al.l 
l200lf ) for the dark matter halos, and the lin e ar an d non- 
linear bias from the work of lScoccimarro et al.l (|200ll ') . Since 
our galaxies are a subsample of the total population, with a 
particular colour and luminosity selection, we do not expect 
that they should follow exactly one track of evolution, but 
in any case, our measured bias evolution seems to agree with 
those galaxies living in halos with masses ~ 10- h-^MQ. 

In our measurements and in the halo model, it is ex- 
pected that when a galaxy tracer has a linear bias ~ 1, 
it should have a small but significantly negative non-linear 
bias 62 • In the right panel of Figure[T5]we show the evolution 
of the WiggleZ galaxies' non-linear bias. These have nega- 
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Figure 13. Left: Constraints on tlie bias parameters bi, 62 and erg for tlie 2(,ft = 0.55 WiggleZ redshift slice. Tiie contours represent l-cr, 
2-cr and S-ct joint confidence regions for a two-parameter fit. Rigfit: Dark Matter reduced 3PCF (black thick line), the WiggleZ Q{d) for 
the Z(,ff = 0.55 slice and the biased dark matter Q{d) (dashed line) using the best-fit parameters found in this analysis. 



N 

1.2 



0.6 
0.4, 



1x 10 



13- 




1x 10' 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 o.s 
redshift 

















_1_x10" 














1x10^"2-.-- 
















1 


t 










1x 10^^ 






1 
































i 


i. 

















0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

redshift 



Figure 15. Evolution of the bias parameters. Left: Evolution of linear bias: triangles are best fit parameters from our WiggleZ regions. 
Lines are halo model prediction of bias for halos of masses A/^ = 10^^ h~^MQ (green dashed-dotted line), = 10^^ h~^MQ (blue 
solid line) and M/^ = 10^^ h~^MQ (red dashed line). Right: evolution of the non-linear bias parameters. 



tive values, and their trend agrees with what is expected of 
~ 10^^ /i-^Mq halos. A more detailed analysis using Halo 
Occupation Distribution models is needed to have a com- 
plete picture of how WiggleZ galaxies populate dark matter 
halos; this is beyond the scope of this paper. 



5.3.2 Evolution of cosmic growth 

In Figure [in] we plot our measurements of as as a function 
of redshift from the WiggleZ survey data. In linear theory, 
the value of erg is calculated as 

a^{R = 8, z) = I ^ \W{k, R = 8)f Puuik, z) (20) 

where W^(fc,i?) is the Fourier transform of a tophat win- 
dow of radius R = 8 h^^Mpc. The linear power spectrum 
evolves as Piin{k,z) oc [D{z) / D{z*)f' P(k, z*) , where z, is a 
reference redshift (e.g., the redshift of recombination) and 
D[z) is the linear growth factor, obtained from the solution 



to the linearized equations of motion of prim ordial overden- 
sities (|Peebleslll980l : lBernardeau et al.ll2002l ). The evolution 
of the linear growth factor depends on th e parameters of the 
cosmological model jLahav fc Sutdl2004l ). Thus, we obtain 

as{z) = ^^^as{z = 0). (21) 

Therefore, asi^z) measurements from the 2PCF and 3PCF 
can be used to study the evolution of the linear growth fac- 
tor. 

As predicted by the standard cosmological model, the 
value of crg(«) we measure decreases at earlier times, in 
agreement with the WMAP5 cosmological parameters. As- 
suming a flat ACDM model with Q,m = 0.27 we find that 
when extrapolated to the present epoch, erg (2 = 0) = 
0.791 

ao7- O^'^ results also agree with the latest estima- 
tion of (jg(z) using B OSS/SDSS-LRGs passive galaxies from 
eiro et aD ()2012l ): modelling the evolution of the linear 
bias for their particular population they find similar values 
to ours. However, in our work we need to make no assump- 
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Figure 16. Evolution of crs{z). Blue symbols correspond to the 
estimates of o-8(2cfl) from the 2PCF and 3PCF of WiggleZ galax- 
ies at different redshifts, marginalizing over the linear and non- 
linear bias parameters. The blue solid line corresponds to the evo- 
lution of crs(z) in a flat WMAP5-ACDM universe, with = 0.27 
and crg{z = 0) = 0.812; the blue shaded region corresponds to 
combined WMAP5 errors. The red dashed line corresponds to 
the evolution of erg (2) in a flat Qrn = 1 CDM universe. The black 
dash-dotted line shoves the evolution of a flat fl\ = 1 universe. 
All models are normalized at the epoch of recombination. 



tions about the evolution of the bias, just in the validity 
(range of scales) of the empirical bias model we adopt. We 
also find agreement with other mea surements of as from the 
3PCF of the SDSS LRG sample jMarfnl I2OI1I ): they find 
as{z = 0.35) = O.BSIq qj, consistent with our measurements 
at the same effective redshift. 

We also plot in Figure [T6]the evolution of as in two dif- 
ferent models of fiat universes, an Einstein-de Sitter model 
(fiat Qrn = i) and another one with n o matter content, just 
cosmological constant. As WMAP5 (jKomatsu et al.ll2009l ) 
measured the amplitude of CMB perturbations at the epoch 
of recombination z ~ 1100, we normalize as at that red- 
shift for the three cosmologies shown here. We find that an 
Einstein-de Sitter universe is disfavoured (Ax^ = 3.5 for one 
parameter fit) when combining measurements of the three 
WiggleZ redshift samples, as well as spacially flat models 
with low matter content (Ax^ ~ 45.01 for Q,m = 0.01). 



clustering, reflecting the morphology of the clustering large- 
scale structures, i.e. the 'cosmic web'. 

• Using a simple local bias parameterization along with 
an empirical treatment of redshift-space distortions of the 
correlation functions, we get constraints on the bias param- 
eters as a function of redshift. Our estimation of the lin- 
ear bias agrees with evolution of dark matter halos ~ 10^^ 
h-'^MQ. 

• For all our redshift samples, we detect a significant non- 
zero (negative) non-linear bias, also consistent with the mod- 
els for the non-linear bias evolution of halos of masses ~ lO^'^ 
h-^MQ. 

• We also constrain the evolution of as with redshift, and 
by extension, the evolution of the linear growth factor. We 
find that our measurements are consistent with the predic- 
tions of a WMAP5 ACDM concordance cosmology and with 
measurements from other methods and observables. 

The improvement in the measurements of the higher- 
order correlations in the last ten years has been dra- 
matic, and it is remarkable that we can now measure the 
3PCF using galaxies up to redshift 2 ~ 1. Although the 
signal-to-noise ratio of the WiggleZ 3PCF is weaker than 
lower-redshift measurements from the SD SS Main Sample 
(iMcBride et al.ll201ll: iMcBride et al.|l201ll ) and SDSS LRG 
sample (|Kulkarni et al.l 1200/1 : iMarfnl 120111 ) . we nonetheless 
have extended the utility of higher-order correlations func- 
tions to 2 1, using the WiggleZ survey data. We note 
that using jack-knife resampling for error esti mation proba - 
bly overestimates the variance on large scales (|Man'nll201ll ): 
our measurements would be improved by the availability of 
mock galaxy catalogues. At the same time, with improved 
statistics we need to improve our modelling of redshift space 
distortions and small-scale effects in order to extract as much 
information as possible from the higher-order correlations. 

Also, improved modelling would help to combine 3PCF 
measurements wi th other observables suc h as clustering 
and lensing (see iMandelbaum et al.l 120121 ') , 2-dimensional 
RSD 2-point stati stics, or the passive galaxies method 
(|Toieiro et al.ll2012l V 

In the near future, with improved measurement tech- 
niques and with bigger surveys, we will be able to use these 
techniques to measure the growth factor accurately and dis- 
cri minate between ACDM model and modified gravity mod- 
els (|Linder fc Cahnll200'iD . 



6 SUMMARY AND CONCLUSIONS 

We have measured the redshift space two- and three- point 
clustering statistics for the WiggleZ galaxies and obtained 
constraints on the linear and non-linear galaxy bias and the 
cosmological parameter ag at three different epochs. Our 
results can be summarized as follows: 

• We obtain significant measurements of the WiggleZ 
3PCF, recovering its shape dependence on large scales, span- 
ning a wide redshift range for all regions and subregions (in 
angle and redshift) of the galaxy sample. 

• These measurements are in agreement with standard 
models for the growth of structure driven by gravitational 
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APPENDIX A: WIGGLEZ 3PCF 
MEASUREMENTS OF ISOSCELES AND VERY 
ELONGATED TRIANGLES 

In Figure lAll we plot results of the WiggleZ correlation 
functions measured for isosceles configurations, where u = 
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Figure Al. The WiggleZ combined correlations of isosceles triangles (n = 1) in redshift-space. We show the absolute 2PCF |5(s3(6))| 
(left), absolute connected 3PCF n = l,0)\ (middle), and reduced 3PCF Q{s,u = I, 6) (right) of WiggleZ galaxies in the z^g = 0.55 
slice (blue triangles), in the slice at z^g = 0.35 (green dashed line) and in the Zgg = 0.68 slice (red dashed-dotted line). Errors have been 
determined by jack-knife resampling. 



S2/S1 = 1, for the three redshift slices we studied. In these 
configurations, the third side length runs from S3 ~ on 
very small angles ^ ~ 0, where the bias model we use is no 
longer valid due to high nonlinearities. Since the third tri- 
angle side S3 covers a large range of scales, for purposes of 
plotting we show the absolute values of the 2PCF and con- 
nected 3PCF on a logarithmic scale, while Q{9) is shown on 
a linear scale (which can take positive or negative values). 
We can observe first that in general the errors in these mea- 
surements are smaller compared with the ones we showed 
in Figure 5 for the u = 2 configurations. As in the m = 2 
configurations, there is no significant evolution in the am- 
plitude of the correlation values. The 3PCF of equilateral 
triangles (u = 1, 6 ~ 1) is small and even negative, as ex- 
pected when galaxies cluster in filamentary structures on 
the largest scales. 

In Figure IA2I we plot results of the WiggleZ correla- 
tion functions measured on very elongated configurations, 
where u = 3, for the three redshift slices we studied. For 
these configurations the signal-to-noise ratio is much smaller 
than in other u configurations, specially on large 9. For the 
Zcff ~ 0.35 redshift slice (green dashed line), ( and Q are 
poorly measured on the largest scales, justifying our deci- 
sion to use a maximum separation S3, max smaller than that 
adopted for the other two redshift slices. Notice that we 
reach scales where the BAO features could in principle be 
observed S3 ~ 100 /i~^Mpc, but the WiggleZ low galaxy bias 
makes it difficult to achieve a significant detection that could 
be used to constrain co smological parameters (as claimed by 
iGaztanaga et al. I I2OO8I . for SDSS LRGs). 
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Figure A2. The WigglcZ combined correlations of very elongated triangles {u = 3) in redshift-space; for the 2PCF ^(53 (^)) (left), 
connected 3PCF Ci^i'^ — 3)^) (middle), and reduced 3PCF Q{s^u = 3,^) (right) of WiggleZ galaxies in the z^ff = 0.55 slice (blue 
triangles), in the slice at z^ff = 0.35 (green dashed line) and in the ZeS = 0.68 slice (red dashed-dotted line). Errors have been determined 
by jack-knife resampling. 



